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The pressure of suspensions of self-propelled objects is studied theoretically and by 
simulation of spherical active Brownian particles (ABP). We show that for certain 
geometries, the mechanical pressure as force/area of a conhned systems can equally 
be expressed by bulk properties, which implies the existence of an nonequilibrium 
equation of state. Exploiting the virial theorem, we derive expressions for the pressure 
of ABPs conhned by solid walls or exposed to periodic boundary conditions. In 
both cases, the pressure comprises three contributions: the ideal-gas pressure due 
to white-noise random forces, an activity-induce pressure (’’swim pressure”), which 
can be expressed in terms of a product of the bare and a mean effective propulsion 
velocity, and the contribution by interparticle forces. We hnd that the pressure of 
spherical ABPs in conhned systems explicitly depends on the presence of the conhning 
walls and the particle-wall interactions, which has no correspondence in systems with 
periodic boundary conditions. Our simulations of three-dimensional APBs in systems 
with periodic boundary conditions reveal a pressure-concentration dependence that 
becomes increasingly nonmonotonic with increasing activity. Above a critical activity 
and ABP concentration, a phase transition occurs, which is rehected in a rapid and 
steep change of the pressure. We present and discuss the pressure for various activities 
and analyse the contributions of the individual pressure components. 
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I. INTRODUCTION 


Living matter composed of active particles, which convert internal energy into systematic 
translational motion, is a particular class of materials typically far from equilibrium. Ex¬ 
amples range from the microscopic scale of bacterial suspension^ to the macroscopic scale 
of flocks of birds and mammalian herds.*^ Such active systems exhibit remarkable nonequi- 
librium phenomena and emergent behavior like swarming,®^ turbulence,® activity-induced 
clustering and phase transitions,®®' and a shift of the glass-transition temperature.®®' The 
nonequilibrium character of active matter poses particular challenges for a theoretical de¬ 
scription. In particular, a thermodynamic description is missing, but would be desirable in 
order to be able to dehne elementary thermodynamic variables such as temperature or chemi¬ 
cal potential.'^ Considerable progress in the description of activity-induced phase transitions 
has been achieved recently.'^2®H23 

Pressure in active fluids attracted considerable attention lately, because a mechanical 
pressure as force per area can be defined far from equilibrium and might be useful in the 
strive to derive a nonequilibrium equation of state.®®I2S12Hl Interestingly, simulation studies 
reveal such an equation of state, with a van der Waals-type pressure-volume phase diagram 
typical for gas-liquid coexistence.'^^®®®' 

In thermal equilibrium, pressure can equivalently be dehned in various ways. In addition 
to the mechanical dehnition, pressure p follows as derivative with respect to the volume 
V of the Helmholtz free energy F, i.e., p = —{dF/dV)T,Ni at constant temperature T and 
particle number Moreover, the virial theorem can be exploited.'^®^ These two identical 
formulations relate the surface pressure to bulk properties. However, for nonequilibrium 
systems, it is a prior not evident how to relate the surface pressure to bulk properties, as it 
has to be kept in mind that there is usually no free-energy functional. By now, pressure in 
active fluids has been calculated using various virial-type expressions ,'I^®®' typically without 
a fundamental derivation. Basic relations are provided for two-dimensional fluids in Ref. EH 
starting from the mechanical dehnition of pressure and by exploiting the Fokker-Planck 
equation for the phase-space probability distribution. 

Here, we apply the virial theorem to derive expressions for the pressure in active huids. 
Both, huids conhned between solid walls and exposed to periodic boundary conditions are 
considered. As an important hrst step, we demonstrate that the mechanical pressure of a 
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confined fluid can equivalently be represented by the virial of the surface forces for particular 
geometries such as a cuboid and a sphere. This has important implications, since it provides 
a relation of the mechanical pressure with bulk properties of the fluid.^ Based on this 
expression, we derive an internal expression for the pressure, which comprises the ideal gas 
pressure, an swim pressure due to particle propulsion, and a contribution by the interparticle 
interactions. As a result, we hud that the swim pressure of the spherical ABPs explicitly 
depends on the presence of conhning walls and the particle-wall interaction, in contrast to 
results presented in Ref. |2I1 but in agreement with more general considerations of anisotropic 
bodies.l^^® Systems with periodic boundary conditions naturally lack such a term. Hence, 
walls can induce features in conhned active systems, which will not appear in periodic 
systems. 

By Brownian dynamics simulations, we determine pressure-concentration relations for 
various Peclet numbers of three-dimensional periodic systems of active Brownian particles. 
As already previously observedj ^^*^^ * ^^ * ^^ we hnd an increasingly nonmonotonic behavior of 
the pressure with increasing propulsion velocity. Above a certain velocity, we also hnd a 
van der Waals-loop-type relation. Most importantly, however, we hnd an abrupt change 
of the pressure for phase-separating huids. Thereby, the magnitude of the pressure jump 
increases with increasing propulsion velocity. Hence, we argue that the nonmonotonic be¬ 
havior indicates increased density huctuations of the active huid, but a phase transition is 
rehected in a more pronounced abrupt change of the pressure. A detailed analysis of the 
individual contributions to the pressure shows that the main ehect is due to the activity 
of the colloidal particles, which gives rise to a so-called ’’swim pressure” as introduced in 
Ref. |26l However, also the interparticle-force contribution exhibits an abrupt increase above 
the phase-transition concentrations. 

The paper is organized as follow. In Sec. [TT| the model of the active fluid is presented. 
The dynamics of ABPs in a dilute system is analysed in Sec. |Hl|for conhned and unconhned 


particles. Expressions for the virial pressure are derived in Sec. IV for conhned and periodic 
systems. Simulations results for periodic systems are presented in Sec. |V| and our hndings 


are summarized in Sec. VI Details of the stochastic dynamics of the ABPs are discussed in 
App. @ and 1^ A derivation of the equilibrium pressure via the virial theorem is presented 
in App. 0 Appendix establishes the relation between the mechanical pressure and the 
surface virial for active systems. 
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II. MODEL 


We consider a system of N active objects, which are represented by hard-sphere-like 
particles of diameter a. A particle i is propelled with constant velocity no along its orienta¬ 
tion vector e, in three dimensions (3D).Its translational motion is described by the 
Langevin eqntion 

f'i(t) = voBiit) + - {Fi{t) + Ti{t )), (1) 

7 

where r, is the particle position, r, the velocity, Fi the total force on the particle, 7 the 
translational friction coefficient, and Ej a Ganssian white-noise random force, with the mo¬ 
ments 


(r*(t)) = 0, 


{Ta,iit)T= 2'ykBT6c,p6ij6{t - t'). 


( 2 ) 

( 3 ) 


Here, T denotes the temperatnre in eqnilibrinm, ks the Boltzmann constant, and a,(3 ^ 
{x,y,x}. The the frication coefficient 7 is related to translational diffnsion coefficient 
via = ksT/'-j. The interparticle force is described by the repnlsive (shifted) Yukawa-like 
potential 


U{r) 


' 

00 , r < a 


ksT < 


u{r) 


V 


— u{rc), r > a ■: 
0 , r > Tc 


( 4 ) 


where u{r) = 5/e j{r — a). Here, I is the interaction range and rc is the upper cut-off 

radius. The orientation performs a random walk according to 


ei{t) = ei{t) X (5) 

where r 7 j(t) is a Gaussian white-noise random vector, with the moments 

{V^) = 0 , ( 6 ) 

= 2Dr5af)5ij5{t - t'). (7) 

Equation ([^ is a stochastic equation with multiplicative noise. An equivalent representation 
for polar coordiantes with additive noise is outlined in App. The translational and rota¬ 
tional diffusion [Dr] coefficients of a sphere are related via The importance 
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of noise on the active-particle motion is characterized by the Peclet number 


Pe = 


vq 

aDr 


( 8 ) 


We study systems in the range of 9 < Pe < 300. Up to = 1.5 x 10^ particles are 
simulated in a cubic box of length L. The density is measured in terms of the global packing 
fraction 0 = Tia^N/ (6U), where V = is the volume. A natural time scale is the rotational 
relaxation time = l/{2Dr). For the potential Q, we use the parameters a/I = 60 and 
Tr = 1.083(7. 


III. DYNAMICS OF ACTIVE BROWNIAN PARTICLES 

A. Unconfined Particles at Infinite Dilntion — Diffnsion 

The equations of motion for a single, force-free particle in an inhnite or periodic system 
follows from Eq. ([^ with Fi = 0. Multiplication of the remaining equation by ri(t) leads to 

^ ^ 

where, (...) denotes the ensemble average. In the asymptotic limit t —)■ oo, the average on 
the left-hand side is the mean square displacement 

(n(i)2) = (r.(0)"> + 6Bp (10) 

for the homogeneous and isotropic system. To evaluate the terms on the right-hand side, 
we use the formal solution 

ri{t) = ri(-oo) + J (voei{t') + dt' ( 11 ) 

of Eq. ([^. Since rji and Pj are independent stochastic processes, the terms on the right-hand 
side of Eq. ([^ reduce to 

{ei{t) ■ ri{t)) = Vo j {ei{t)-ei{t'))dt'= (12) 

1 /■* 

(R(t) ■ r,{t)) = - (R(t) ■ R(F)) dt’ = SksT. (13) 

T J —OO 

Here, we use the correlation function (cf. App. 0 

(e,(^)■e,(0))=e-2^’'^ (14) 
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and the definition 


1, t G (— 00 , t) 



1 

= 12 - * = ‘ 


(15) 


— 00 


0 , t G (t, 00 ) 


Hence, Eq. ([^ yields the diffnsion coefficient 



(16) 


in the asymptotic limit t —)• 00 . This is the well-know relation for the mean sqnare 


displacement of an active Brownian particle in three dimensions.^ However, we have de¬ 
rived the expression exploiting the virial theorem.^ More general, in d dimensions follows 
Dl^ = + v^/d^d — The expression has been confirmed experimentally for 2D 


synthetic microswimmers in Refs. [33 and 133 

B. Confined Particles — Mean Particle Velocity 

For ABPs confined in a volnme with impenetrable walls, the force on particle i is 

N 

F.(i)= + (17) 

i=i 

where the F^j = Fij{ri — Vj) are pairwise interparticle forces and Ff is the short-range 
repnlsive force with the wall. We assnme that the wall force points along the local surface 
normal n, i.e., Ff = — F/rij, where nj = n(rj) points outward of the volume. The prime at 
the sum of Eq. (0 indicates that the index j = i is excluded. 

By multiplying Eq. ([^ with the orientation vector ei{t) and averaging over the random 
forces, we obtain 



(18) 


The expression 



( 19 ) 


i=l 


6 





can be considered as the average particle velocity along its propnlsion direction in the inter¬ 
acting system.l^ Hence, with Eqs. (17) and (18), the propnlsion velocity can be expressed 


as 

N N 1 ^ 

^ = >'» + Vv ^ ^ ^ ■ 

' i=i j=i ' i=i 

Aside from the wall term, a similar expression has been derived in Ref. 
interaction term can be written as 


N N 


N N 


E' § E E' ■ («• - ■ 


( 20 ) 

The pairwise 

( 21 ) 


i=l j=l 


i=l j=l 


As a conseqnence, we obtain a contribntion to the velocity only when jj' ej. 

The velocity perpendicular to a wall of a particle interacting with a wall is zero. Hence, 
the force (Ff) of such a particle, averaged over the random force, is 


{F‘) = 7tio (e, ■ Tii) 


( 22 ) 


for a dilute system. Thus, the wall force yields a contribution to the mean propulsion velocity 
as along as / n^. In case of a dilute system, the mean velocity reduces then to 

.. N N 

^ (1 - • ( 23 ) 

^ i=i i=i 

Hence, the presence of walls reduces the mean propulsion velocity. This is in accord with 
the considerations in Ref. |25] for anisotropic particles, but an extension of the calculations 
of Ref. [2T]for spherical ABPs. Naturally, the velocity is zero for particles adsorbed at a wall 
and an orientation vector parallel to n. 

Generally, the average propulsion velocity of a particle, which is part time in the bulk, 
with probability Ptuik, and part time at the wall, with probability 1 — Ptuik, is 

Vi = voPbuik + no (l - ((e* ■ riif)) (1 - Pbuik)- (24) 


Thus, the distribution function of the angle between e* and n* is required, as well as the 
probability distribution to hnd a particle in the bulk to determine the average propulsion 


velocity. Naturally, the average in Eq. (23) can be calculated by the distribution function 
of the particle position and the propulsion orientation e^, which follows from the respective 
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Fokker-Planck equation with appropriate boundary conditions 


More general, multiplication of Eq. ([^ by the velocity hj(0) yields 

{ei{t) ■ ri{0)) = (ei(0) ■ ri(0)) (25) 

within the Ito interpretation of the stochastic integral.®® Then, averaging over all particles 
gives 

1 ^ 

( 26 ) 

i=l 

This relation describes the average decay of the correlation between the initial velocities 
rj(0) and the time-dependent orientations ej(f). Not surprisingly, the correlation function 
decays exponentially, because the dynamics of the orientation ei{t) is independent of the 
particle position and any particle interaction. 

IV. VIRIAL FORMULATION OF PRESSURE 

The expression of the pressure depends on the boundary conditions. Here, we will discuss 
a system of ABPs conhned in a cubic simulation box with either impenetrable walls or 
periodic boundary conditions. 

A. Confined System 

Multiplication of the equations of motion ([^, with the forces 10. by ri{t) yields 

quo {ei{t) ■ n(f)) + (ri(f) ■ ri{t)) + {Fi ■ n) = 0, (27) 

since the mean square displacement is limited for the conhned system, and hence, 

its derivative vanishes in the stationary state. The force contribution can be written as 

N N N N 

E ■’■*> = 2 E E' <-^6 ■ (n - rj)) + E {F‘ ■ n) . ( 28 ) 

2=1 2 j = l 2=1 

The hrst term on the right-hand side is the virial of the internal forces and the second 
term that of the wall forces. The latter needs to be related to the pressure of the system. 




which is the major step in linking the mechanical (wall) pressnre to bnlk properties. The 
traditional derivation for equilibrium systems is briefly outlined in App. 0 This derivation 
assumes a homogeneous pressure throughout the volume. In contrast, active systems can 
exhibit pressure inhomogeneities depending on the shape of the confining wall,*^^*^ and thus, 
the derivation does not apply in general. So far, there is no general derivation or proof of 
the equivalence of the external virial and the pressure; there may be no such universal 
equivalence, since there is no thermodynamic potential for an active system. 

However, for special geometries, we can establish a relation between the wall forces and 
pressure. This particularly applies to cuboidal and spherical volumes. As shown in App. 
for such geometries the average wall pressure is given by 


N 


3p‘V = -J2 (U* ■ n ), (29) 

i=l 

where we introduce the external pressure to indicate its origin by the external forces. As 


an important result, summation of Eq. (27) over all particles yields the relation between 
pressure and bulk properties 

N ^ N N ^ 

3p^V = SNksT + 7^0 (e* + ■ i^i - rj )). (30) 

2 = 1 2=1 j = l 


Naturally, p* = p®. However, we introduce the internal pressure p* to stress the different 
nature in the pressure calculation by bulk properties. The first term on the right-hand 


side of Eq. (30) is the ideal-gas contribution to the pressure originating from the random 


stochastic forces Fj (Eq. (13)). Without such forces, our system reduces to that considered 
in Ref. [T31 The other two terms are the virial contributions by the active and interparticle 
forces. Following the nomenclature of Ref. ESI we denote second term on the right-hand side 


as ’’swim virial”. It is a single-particle self contribution to the pressure, in contrast to the 
virial contribution by the interparticle forces.^ Hence, the swim virial resembles similarities 
with a contribution by an external field rather than interparticle forces. 

The contributions to the swim virial can be expressed as 

nt poo 


{ei{t) = / {ei{t) ■ dt' = / (ei(f') ■ ri(0)) dt' 


(31) 


in the stationary state and with (rj(—cx)) ■ ej(f)) = 0. Together with Eq. (26), Eq. (30) can 
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be expressed as 


AT . N N 

3p-V^3NksT+2^v+lY.T!(^‘r 


2Dr 


rj)) 


(32) 


i=l j=l 


in terms of the mean propulsion velocity v (Eq. 20). A similar expression is obtained in 
Ref. EH but in a very different way. However, in contrast to Ref. EH our mean velocity 


(Eq. 20) comprises interparticle as well as wall contributions. Hence, the wall interactions 
are not only manifested in the pressure, but appear additionally in correlation functions, 
i.e., (ej(f) ■ Viit)). This a a major difference to passive systems, were Vq and hence v are 
zero. 


In the case of zero interparticle forces, the pressure Eq. (32) reduces to 

N 


Sp^V = SNkeT 1 + 




0 


6DrD^ 


■ 


(33) 


2 = 1 


The hrst term on the right-hand side corresponds to the ideal (bulk) pressure 


pi'i = 


NksT 

V 


1 + 




(34) 


of an active system.^^IESElllI] second term accounts for wall interactions, and is respon¬ 
sible for various phenomena such as wall accumulation.^I^IiSEZllI] 


B. Periodic Boundary Conditions 

In systems with periodic boundary conditions, there is no wall force but in addition 
to interactions between the N ’’real” particles, interactions with periodic images located at 
Vi + Rn appear. For a cubic simulation box, the vector R^ is given by 

R^ = (35) 

with ria G Z. The generalization to a rectangular noncubic box is straightforward. For a 
pair-wise potential Uij, the potential energy is then given by 

i=l j=l n 

In general, this involves inhnitely many interactions.®^ Here, we will assume short-range 
interactions only. This limits the values Ua to Uq ± 1 for a particle in the primary box, which 
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corresponds to an application of the minimum-image convention.®! The forces in Eq. Q 
then given by 


are 


N 


Fi — Fij [vi Vj Rr, 


(37) 


j=l n 


Multiplication of Eq. ([^ by and summation over all particles leads to the relation 


^^0 


iV iV -‘v 


Z=1 


27 

1 


(38) 


i=l j=l n 

N 1 ^ ^ 

2 = 1 


7t=i 


with the abbreviation Fij{ri — Vj — R^) = F^. By following the trajectories of the particles 
through the inhnite system, i.e., by not switching to a periodic image in the primary box, 


when a particle leaves that box, {riity) is equal to the mean square displacement of Eq. (10), 
but with a translational diffusion coefficient Dt of the interacting ABP system. Moreover, 


Eq. (13) applies and, hence. 


N 


N N ^ 

3'yNDt = SNksT + 7^0 5^ (e* + - ^^■)> • ^39) 


2 = 1 


2=1 j=l n 


In the force-free case, i.e., F^ = 0, this relation is identical with Eq. (|16|). More general. 


Eq. (39) provides the self-diffusion coefficient of active particles in a periodic system inter¬ 


acting by pair-wise forces. Exploiting Eqs. (26) and (31), the active virial can be expressed 


by the mean propulsion velocity v, and the diffusion coefficient becomes 

N N 

' 2=1 j=l n 

with the mean propulsion velocity in a periodic system 

N N 


(40) 


(41) 


2=1 j=l n 


After insertion of the velocity n, Eq. (40) provides three contributions to the diffusion coeffi¬ 


cient: (i) the diffusion coefficient, Eq. (16), of non-interacting active Brownian particles, (ii) 
a contribution due to correlations between the forces F^ and the orientations e*, and (iii) 
a contribution from the virial of the inter-particle forces. This clearly reflects the explicit 
dependence of the diffusion coefficient on the orientational dynamics of the particles. 
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In contrast to confined systems, the multiplication of the equations of motion by the 
particle positions does not directly provide an expression for the pressure of a periodic 
system. Here, additional steps are necessary. To derive an expression for the pressure, we 
subtract the term 


N N 




(42) 


i=l j=l n 


from both sides of Eq. (39). Then, we are able to dehne an external {p'^) and internal (p*) 


pressure of the periodic system according to 

D 1 ^ ^ 

3pV = SNk^T- - - 5^ 5^' 5^ (jpy ii..) . 

* i=l j=l n 


(43) 


N ^ N N ^ 

3 p^V = SNksT + 7^0 5 ^ (e* + -Rrr)), ( 44 ) 


2=1 


2=1 j=l n 


as for passive systems in Refs. [301 - I321 Naturally, both expressions yield the same pressure 
value, i.e., p® = ph 

For notational convenience, we write the expressions for the pressure as 


p*^ = p'^ + p®'^, (45) 

p* = p° + p"^ + p*'^, (46) 


with the abbreviations for the individual terms of Eqs. (43) and (44): 

NkeT 


p° — 

^ AT 

P =777 2^ (e* ■ rt) = ^T^yr^VoV, 


3V 

2=1 

, NksT Dt 


6VDr 


p 


p^f = 


V DT 

N N 


6V 


i^n)) , 


1=1 j=l n 


N N 

7 = --EE'ET7--fu). 


6V 


i=l j=l n 


(47) 

(48) 

(49) 

(50) 

(51) 


The expressions Eqs. (43) and (44) are extensions of those of a conservative passive sys¬ 
tem, where Dt = 0 and Vq = 0, with a similar dehnition of external and internal pressure.®^!] 
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For both cases, the dehnitions appear ad hoc. Note that for a passive system typically New¬ 
tons’ equations of motion are considered. However, for the passive system, the expression for 
the internal pressure can be derived from the free energy,®! as illustrated in the appendix of 
Ref. HU which underlines the usefulness and correctness of the applied procedure for passive 


systems. The similarity between the internal pressure expressions Eqs. (30) and (44) for 


the conhned and periodic systems, respectively, supports the adequateness of the applied 
procedure in the pressure calculation. 

The individual terms can be understood as follows. 

External pressure — The pressure contribution accounts formally for the forces 


across faces of the periodic system, similar to the expression (29) of a conhned system 


There is only a contribution for non-zero vectors Rn, i.e., for interactions between real and 
image particles.®* The diffusive dynamics yields the pressure contribution p‘^ ~ Dt/D^. This 
term is not present in a conservative passive system and is a consequence of the stochastic 
forces Fj .ESESI por passive, non-interacting particles, i.e., Uq = 0 and U of Eq. (36) is zero. 


Dt = and the pressure reduces to the ideal gas expression p® = p^. For non-interacting 
but active systems, Dt is given by Eq. (16) and p® = p®'’* (Eq. (34)). 


Internal pressure p® — The internal pressure comprises the ideal gas contribution p°, the 
swim pressure p^ due to activity, and the contribution p®'f by all the pair-wise interactions in 
the system within the minimum-image convention. As for the virial of the conhned system 


(32), the swim pressure p® Eq. (48) can be expressed by the mean propulsion velocity v 


(Eq. (41)). The forces can be strongly correlated with the propagation direction Cj, 
specihcally when the system phase separates as in Ref. dni For non-interacting particles or 


at very low concentrations, the contribution of the intermolecular forces vanishes and the 


pressure is given by Eq. (34). 


There are various similarities between the expressions for the pressure of a conhned and a 


periodic system. Formally, Eqs. (32) and (44) are similar. However, specihcally the meaning 


of the mean propulsion velocity v is very diherent. The presence of the wall term in Eq. (20) 
can lead to wall-induced phenomena, which are not present in periodic systems. Hence, 
results of simulations of periodic and conhned active systems are not necessarily identical. 
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FIG. 1. Schematic phase diagram of active Brownian particles in terms of the Peclet number Pe 
and the packing fraction (j). Symbols correspond to the Peclet numbers for which the pressure is 
calculated. The shaded areas indicate the liquid-gas (light green) and the crystal-gas (light yellow) 
coexistence regimes. Further details can be found in Ref. [191 

V. SIMULATION RESULTS 

Simulations of active Brownian spheres with periodic boundary conditions (cf. Sec. 0 
reveal an intriguing phase behavior and a highly collective dynamics, as discussed in Ref. [19] 
for three dimensional systems. At low concentrations, the ABPs exhibit a uniform gas¬ 
like phase. However, above a critical, Peclet-number-dependent concentration, the isotropic 
system phase separates into a dilute, gas-like phase, and a dense, liquid-like phase.^^^'^ 
Figure [T] displays the phase diagram for the considered concentrations and Peclet numbers. 
Above Pe ~ 30, we hnd a wide range of concentrations, over which the system phase 
separates. 
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FIG. 2. Translational diffusion coefficients of ABPs for the Peclet numbers Pe = 9.8 (purple), 29.5 
(olive), 44.3 (green), 59.0 (blue), and 295.0 (red). The values are normalized by the respective 
diffusion coefficients ( |16[ ) of an individual ABP. The inset shows the same data on linear scales. 
The solid line (black) indicates the dependence Dt/Dl^^ = 1 — 2(j). 


As is evident from Eqs. (43) and (44), diffusion contributes to the pressure. Hence, we 
address the dynamics of the ABPs in terms of the diffusion coefficient in the following. 
Subsequently, we discuss the dependence of the pressure and its individual contributions on 
Peclet number and concentration. 


A. Diffusion coefficient of ABPs 

We determine the diffusion coefficient of the ABPs by the Einstein relation, i.e., we 
calculated the particle mean square displacement and extracted from the linear long-time 
behavior the translational coefficient Dt- Figure displays the obtained values as function 
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of the packing fraction for various Peclet numbers. For small packing fractions (0 < 0.3), 
Dt/Dl'^ decreases linearly with increasing concentration. The decay is reasonably well de¬ 
scribed by the theoretically expected relation 1 — 20 for colloids in solution.®^ However, a 
better agreement is obtained for the relation 1 — 2.20, i.e., a 10% large factor. Remarkably, 
the decrease of Dt for 0 < 0.3 is essentially independent of the activity. Signihcant devia¬ 
tions between the individual curves are visible for higher concentrations and larger Peclet 
numbers, where a phase separation appears. At the critical concentration, the diffusion 
coefficients for Pe > 30 exhibit a rapid change to smaller values. A further increase in con¬ 
centration implies a further slowdown of the dynamics. Interestingly, at high concentrations, 
the larger Peclet number ABPs diffuse faster than those with a small Pe. We attribute this 
observation to collective effects of the ffuids exhibiting swirl- and jet-like structures.^ An¬ 
other reason might by an activity induced fluidization, since the glass transition is shifted 
to larger packing fractions.!^ 


B. Pressure of ABPs 


Figure [^displays the density dependence of the total pressure p = for the consid¬ 


ered Peclet numbers. We normalize the pressure by the ideal gas value (34) and multiply by 


the packing fraction 0 in order to retain the linear density dependence. As expected, both. 


Eq. (45) and Eq. (46) yield the same pressure values. We observe a pronounced dependence 


of the pressure on the Peclet number. Note that the pressure for the largest Peclet num¬ 
ber is highest, because ~ Pe^ for Pe 3> 1. For small Pe, the pressure increases with 
increasing 0. However, already for Pe = 29.5, we obtain a strongly nonmonotonic concen¬ 
tration dependence, where the pressure decreases again for 0 > 0.35. A similar behavior 
has already been obtained in Refs. [13], [ST] and |26l In Ref. |26l the decreasing pressure at 
high concentrations is attributed to cluster formation of the ABPs. For ABPs with Peclet 
numbers well in the phase separating region, we hud a jump of the pressure at the critical 
packing fraction. This behavior has not been observed before, most likely because smaller 
Peclet numbers have been considered only. Such a pressure jump seems not to be present 
for Pe = 29.5, but the Peclet number is close to the binodal (cf. Fig. and a clear phase 
separation is difficult to detect in simulations. 

Interestingly, all the pressure curves exhibit a remarkably similar behavior for Pe ^ 29.5 
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FIG. 3. Pressure in systems of ABPs for the Peclet numbers Pe = 9.8 (purple), 29.5 (olive), 
44.3 (green), 59.0 (blue), and 295.0 (red). The values are normalized by the pressure (34) of 
individual ABPs at the same Peclet number. The black line indicates the dependence = 

(/)(1 — 1.2(/>), and the purple line the relation = (f){l — 0.^ 


and 0 < 0.35 and can very well be fitted by the density dependence 

p = 0(1 - «0) (52) 

in this regime, with k ~ 1.2, as indicated in Fig. For smaller Pe = 9.8, we hnd k ~ 0.8. 
Hence, in any case the pressnre increases linearly for 0 <C 1, as expected, and exhibits 
qnadratic corrections for larger concentrations, which depend on the swimming velocity or 
Peclet nnmber. In Ref. UHl the noneqnilibrinm virial eqnation of state 

p = [1 - (1 - 3/2Pe)0] (53) 

is proposed, adopted to onr dehnition of the Peclet nnmber, and = NksTil — Pe^/2) /V. 
This relation applies only within a certain range of Pe valnes. It certainly acconnts for the 
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FIG. 4. Active pressure Eq. (48) in systems of ABPs for the Peclet numbers Pe = 9.8 (purple), 
29.5 (olive), 44.3 (green), 59.0 (blue), and 295.0 (red). The values are normalized by the pressure 


(34) of individual ABPs at the same Peclet number. 


increase of the prefactor of 0 with increasing Pe, but the asymptotic value for Pe i:§> 1 
does not agree with our simulation result. However, our simulation results for Pe = 9.8 are 


reasonably well described by Eq. (53) at small concentrations. 


The nonlinear concentration dependence of Eq. (52) with a positive k yields a negative 


second-virial coefficient B 2 = —7ra^K/6. In classical thermodynamics, a negative B 2 suggest 
the possibility of a gas-liquid phase transition.^ For the ABPs, the negative B 2 definitely 
indicates large concentration fluctuations, which increase with increasing Peclet number. 
However, in the light of our simulation results, it seems that a phase transition is also 
associated with a pressure jump at the critical concentration. This is more clearly expressed 
by the swim pressure displayed in Fig. Only for Pe > 29.5, we hnd a pressure jump 
and a phase transition. 
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FIG. 5. Pressure and its individual contributions of ABPs for the Peclet number (a) Pe = 9.8 
and (b) Pe = 295.0. The values are normalized by the pressure (34) of individual ABPs at 


the same Peclet number. The green symbols indicate the total pressure Eqs. (45) and (46). The 
red symbols correspond to internal pressure coii|)onents (■) and p*'^ (•) of Eqs. (48) and (50), 


respectively. The blue symbols indicate the pressure components p^ (■) and (•) of Eqs. (49) 
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FIG. 6. Contribution of the interparticle force (Eq. (50)) to the internal pressure of ABPs for the 
Peclet numbers Pe = 9.8 (purple), 29.5 (olive), 44.3 (green), 59.0 (blue), and 295.0 (red). The 
solid line (black) indicates the dependence ~ 


Figure illustrates the contributions of the individual terms of Eqs. (48) - (51) to the 
total pressure for the Peclet numbers Pe = 9.8 and Pe = 295. At small Pe = 9.8 (see 
Fig. I^a)), the terms p® and decrease monotonically with increasing concentration and 
approach zero at large (j). For p^, this is easily understood by the slowdown of the diffusive 
dynamics due to an increasing frequency of particle encounters with increasing density. The 
pressure p'^ is given by 


N N ^ 

/ = ^ ^ . (r, - r,)> (54) 

i=l j=l n 

in terms of the ideal gas and the swim pressure. As for a passive system, the last term on the 
right-hand side of this equation is smaller than zero, and hence, p® > p'^ — p°. The equation 
also shows that the swim pressure is generally different from p'^ and hence is not simply 
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FIG. 7. Contribution of the interparticle force (Eq. (|51[)) to the external pressure of ABPs for the 
Peclet numbers Pe = 9.8 (purple), 29.5 (olive), 44.3 (green), 59.0 (blue), and 295.0 (red). The 
solid lines indicates the dependencies jp^'^ ps 1.35(?!) (purple) and p^f /p^'^ ps — (j)) (black). 

related to diffusive behavior of the ABPs. As far as the force contributions and are 
concerned, they monotonically increase with 0 for Pe = 9.8, until the ABPs undergo a phase 
transition and form a crystalline structure. The gradual increase of pairwise interactions 
with increasing concentration leads to a systematic increase of this pressure contribution. 
The ratio decreases initial {(f) < 0.4) due to the slowdown of the diffusive dynamics, 

but at large (f>, the contributions p®-^ and p*-^ dominate and the respective pressure increases 
with concentration similar to that of a passive system. For Pe = 9.8, p is essentially equal 
to p®-^ for (f> > 0.5. 

The relevance of the individual contributions changes with increasing Peclet number. For 
Pe = 295 (see Fig. [^b)), p'^ and p^ still decrease monotonically with increasing (f>. However, 
the total pressure is essentially identical with p^, i.e., p is dominated and determined by 
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the swim pressure. The contribution pT |;o interparticle interactions is very small and 
is only relevant at high concentrations. The pressure increase for small concentrations, 
passes through a maximum and decreases again for large 0. This nonmonotonic behavior is 
very different from that at low-Peclet numbers. The appearing jumps in p® and p"^ reflect the 
appearing phase transition at a critical density Note that the critical density depends 
on the Peclet number. 

Aside from absolute values, the density dependence of the various pressure components is 
rather similar for densities below the critical density (/>c ~ 0.35. In any case, the pressure is 
dominated by the active contribution as expressed and reflected in the diffusion coefficient. 
Above (pci however, the contribution by the inter-particle forces dominates for Pe = 9.8, 
whereas p is still determined by p® for Pe = 295. Here, the contribution of p^ exceeds that 
of p^'f by far, although the latter increases also with the Peclet number as shown in Fig. but 
p^ ~ Pe^ at small 0 and outweighes the growth of p*-^. Two effects contribute to the decrease 
of the overall pressure. The diffusion coefficient decreases with increasing concentration due 
to an increasing number of particle-particle encounters. However, p"^ decreases faster than 


p^. Via the relation (54), we see that this faster decrease is compensated by an increase in 
the contribution by the interparticle forces. 


Looking at the swim pressure p^ (cf. Eq. (41)), it evidently decreases because the cor¬ 
relations between the forces and the orientations Cj decrease with increasing concentration. 
Here, an important aspect is the appearance of collective effects of the ABPs, specifically 
for (j) > Collective effects lead to a preferential parallel alignment of certain nearby 


ABPs. As indicated in Eq. (21), such an alignment reduces the force-propulsion-direction 
correlations, and hence, the pressure. 

As indicated before, the force contribution to the internal pressure p*-^ depends roughly 
linearly on Pe, as illustrated in Fig. Remarkably, p®-f exhibits an exponential dependence 
on density, namely p*-^ ~ exp (60), over a wide range of concentrations for Peclet numbers 
in the one-phase regime, and up to the critical value 0c for Peclet numbers where a phase 
transition appears. Additionally, the forces exhibit a rapid increase at the critical density. 
We did not analyse in detail whether a jump appears or if pT is smooth and a gradual but 
rapid increase occurs, as might be expected from the data for Pe = 44.3. 

Figure [^displays the contributions of the internal forces (Eq. (51)) to the external pressure 
(Eq. ([45|). The general shape of the curves is comparable to that of the scaled pressure 
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Compared to the latter, exhibits the stronger density dependence ~ 0^ for 
0 1. However, compared to p^'f, the density dependence is weak. Quantitatively, it can 

be described as p®f/p*‘^ = k0(l — 0) with k ~ 1.35 over a wide range of Peclet numbers 
9 < Pe < 30. With increasing Pe, k decreases slightly and we find k ~ 1.1 for Pe = 295. 


VI. SUMMARY AND CONCLUSIONS 

We have presented theoretical and simulation results for the pressure in systems of spher¬ 
ical active Brownian particles in three dimensions. In a first step, we have shown that the 
mechanical pressure as force per area can be expressed via the viral by bulk properties for 
particular geometries, especially a cuboid volume as typically used in computer simulations. 
Subsequently, we have derived expressions for the pressure via the virial theorem of systems 
conhned in cuboidal volumes or with periodic boundary conditions. For the latter, two 
equivalent expressions are derived, denoted as internal and external pressure, respectively. 
As a novel aspect, we considered overdamped equations of motion for active systems in the 
presence of a Gaussian and Markovian white-noise source. The latter gives rise to a thermal 
diffusive dynamics for the translational motion and yields the ideal-gas contribution to the 
pressure. 

The activity of the Brownian particles yields distinct contributions to the respective 
expressions of the pressure. For the external pressure, an effective diffusion coefficient of 
an active particle appears. The internal pressure contains a term, which is related to the 
product of the bare propulsion velocity of an ABP and its effective mean propulsion velocity. 
The latter is determined by the interactions with surrounding ABPs and reduces to the bare 
velocity at inhnite dilution. Similar expressions have been derived in Refs.^^^^^HMI 

Our simulations of ABPs in systems with three-dimensional periodic boundary conditions, 
show that the total pressure increases first with concentration and decreases then again at 
higher concentrations for Peclet numbers exceeding a ceratin value, in agreement with pre¬ 
vious studiesThereby, the pressure exhibits a van der Waals-like pressure-concentration 
curve. This shape clearly reflects collective effects, however, it does not necessarily indicate 
a phase transition. In case of a phase transition, we hnd a rapid and apparently discontinues 
reduction of the pressure at the phase-separation density. 

Quantitatively, the density dependence of the pressure is well described by the relation 
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p ~ 0(1 — K0), where the range over which the formula applies and k depends weakly on 
the Peclet number. For Pe ~ 10, we hnd k ~ 0.8. For Peclet numbers in the vicinity and 
above Pcc, all the considered curves are well described with k, 1.2. Generally speaking, 
our K, values extracted from the simulations are in the vicinity of the value k = 1 proposed 
in Ref. ES 


Considering the various contributions to the pressure individually, we hnd that the in¬ 
ternal pressure is dominated by the contribution from activity for Pe > Pcc, in particular 
for 0 < 0c. For Peclet numbers larger than the critical value, there is signihcant contribu¬ 


tion from the interparticle forces of Eq. (41). This contribution decreases with increasing 
Peclet number, since the swim pressure increase quadratically with the Peclet number and 
the pressure due to interparticle force linearly. However, for moderate Pe, the exponential 
increase of with concentration outweighes the Peclet-number dependence. 


In summary, the pressure of an active system exhibits distinctive differences to a passive 
system during an activity-induced transition from a (dilute) huid to a system with a coexist¬ 
ing gas and a dense liquid phase. While the pressure of a passive system remains unchanged 
during the phase transition, our simulations of an active system yield a rapid decrease of the 
pressure for densities in the vicinity of the phase-transition density. Thereby, the pressure 
change becomes more pronounced with increasing Peclet number. Hence, evaluation of the 
pressure, specihcally the swim pressure, provides valuable insight into the phase behavior 
of active systems. 


We gratefully acknowledge support from the DFG within the priority program SPP 1726 
” Microswimmers”. 


Appendix A: Rotational Equation of Motion 

Equation (j^ is a stochastic equation with multiplicative noise. The Fokker-Planck equa¬ 
tion, which yields the correct equilibrium distribution function, follows with the white noise 
in the Stratonovich calculus.*®®^ 

For the numerical integration of the equations of motion (|^, an Ito representation is 
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more desirable. Introducing spherical coordinates according to 


e = 


cos 93 sin";? 
sin(y 9 sin-d 
cos'd 


\ 

/ 


(Al) 


(here, we consider a particular particle and skip the index i), alternatively the stochastic 
differential equations can be considered 


e = T]^ei) + (A2) 

and yields the same Fokker-Planck equation for the generalized coordinates (p and 
Thereby, the equations of motion of the angles are 


with additive noise, 
and normalization, 
moments 


'd = D cot 'd + , 

1 


(p = 


sin-d 


hi? 


(A3) 

(A4) 


The unit vectors e^, ^ G {(p,'d}, follow by differentiation ~ de/d^ 
The 7]^ are again Gaussian and Markovian random processes with the 


(he) = 0, 

(he(^)he'(^')) = 2T)^%/d(t - t'). 


(AS) 


For simulation, Eq. (A2) has to be discretized. Integration, iteration, and usage of the 


Eqs. (A3) and (A4) yields the difference equation 


e{t + At) = e{t) + e^{t)Ar]^ + e^{t)Ar]^ — 2Dre(t)At. 


(A6) 


The Arj^ are dehned as 


/ t+At 


(A7) 


and are Gaussian random numbers with the hrst two moments 


(Ahe) = 0 

{Ar]^Ar]^i) = 2Dr6^^/At. 


(AS) 
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This is easily show with the help of Eqs. (A5).^As common for the Ito calculus, we replaced 
the values At^I by their averages 20^At in Eq. (A6) 


For practical reasons, it is more convenient to apply the following integration scheme. 
An ’’estimate” e'{t + At) is obtained via 


e'{t + At) = e{t) + e^{t)Ari^ + e^{t)Ar]^. 


(A9) 


This new orientation vector is no longer normalized. Normalization yields e{t + At) = 
e'{t + At)/\e'{t + Af)|, where \e'{t + Af)p = e'{t + At) ■ e'{t + At) = 1 + Arj^ + Arj^. Since 
Arj^ ~ At, we obtain for DrAt 1 

e(t + At) = e(t) + e^{t)Ar]^ + e^{t)Ar]^ 

- \ e{t) (AlO) 

up to order At^/^. Upon replacement of Ar/I by its average, we obtain the desired equation 

«)• 


Appendix B: Orientational Vector Correlation Function 

The autocorrelation function of the orientational vector is conveniently obtained from 
Eq. (|^ for an infinitesimal time interval dt, i.e., the equation 

de = e(t) X dri{t) — 2Dredt (Bl) 

within the Ito calculus. Multiplication with e(0) and averaging yields 

{de{t) ■ e(0)) = ([e(t) x drj{t)] ■ e(0)) — 2Dr {e{t) ■ e(0)) dt. (B2) 

Since e{t) is an nonanticipating function,!^ the average over the stochastic-force term van¬ 
ishes, and we get 

|(e(i)-e(0)) = -2D,(e(«).e(0)), 

which gives 

(e(f) -6(0)) 
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Appendix C: Virial and Pressure in Equilibrium Systems 


To relate the external virial with the pressure for an equilibrium system, we consider a 
surface element AS and the forces exerted on that surface element Ff. The pressure 

p® at the surface element is 


P 


'^s =j2iP‘) = -Y, ■ 


(Cl) 


i&AS ieAS 

On the other hand, the external virial reads as 


ieAS \ ieAS 

Since the surface element is small, the vector is essentially the same for all particles 
interacting with AS and can be replaced by the vector r of that element. The same applies 


to the normal n*, i.e., n* = n, where n is the local normal at AS. With Eq. (C2) and the 
mechanical pressure, we then obtain 


E m = -r nJ2 iF,‘) = -p‘r ■ &S, 


(C3) 


ieAS ieAS 

with AS = nAS. Integration over the whole surface gives 


N 

fp‘r-dS = J2iF‘-n}. 

'' i=l 


(C4) 


Using Gauss theorem, the surface integral turns into a volume integral, i.e.. 


^ r ■ dS = J 'V ■ (pr) d^r = 3p®U. (C5) 

For the last equality, we assume that the pressure is isotropic and homogeneous. Hence ,1^21 

N 

3p‘v = - E {F‘ ■ n ). 


(C6) 


2=1 


Appendix D: Virial and Pressure of Active Particles 
1. Cuboidal Confinement 

We consider particles conhned in a cuboidal volume of dimensions x Ly x L^, with the 
reference frame located in its center (cf. Fig. |^. The total mechanical pressure is given by 
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FIG. 8. Cross section of a cuboid for the calculation of the mechanical pressure. 


the average over the six surfaces, i.e., = {pi + pu + ... + pvi)/Q- For the surface Si, the 

mechanical pressure is 

i 

Since n ■ Vi = for any point on the surface, we introduce unity, i.e., n ■ ri/{Lx/2) = 1 
and obtain 


Pi = 


SrL 


2 2 
E w'■’•<> =E w'. 


(D2) 


with = —Ffn. Hence, the pressure is given by 

N 

3p‘V = - ^ (f;* ■ r ,}, (D3) 

i=l 

as for an equilibrium system. The derivation of this relation is independent of the choice of 
the Cartesian reference frame. At least, as long as the axis unit-vectors are normal to the 
surfaces. 

28 










2. Spherical Confinement 


For a sphere of radius R with the reference frame in its center, the pressure is 


N N 


2=1 


(D4) 


2 = 1 


where S = For every point on the surface applies = i?n, hence, 


N 


= 


SR 




(D5) 


2 = 1 


With SR = 4:71 R^ = 3C, we again hnd the relation (D3). 
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